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This paper is timely for highlighting the benefits of Quasi-Monte Carlo (QMC) in con¬ 
temporary computational statistical methodology. Below we address the question of whether 
there is scope to further reduce the error of QMC estimators. The analysis of QMC used by 
Gerber and Chopin is rooted in the Koksma-Hlawka inequality 
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is a test function of interest, is a point set (or sequence), V{(p) is 


where cp : [0,1]'^ —)■ 

the (Hardy-Krause) total variation and is the (star) discrepancy term that is the 

target of the QMC innovation. Our discussion explores the potential to simultaneously tackle 
the rate constant U(y) in conjunction with the use of QMC methods to tackle 
This direction has received considerably less attention due to typical analytic intractability 
of the rate constant. [Hickernell et al. (2005) showed that classical control variate strategies 
from Monte Carlo (MC) are typically not well-suited to QMC, since the total variation is 
only weakly related to the MC variance that is the target of classical variance reduction 
techniques. Below we hint toward a general strategy to reduce QMC error that targets the 
rate constant directly. 

Following recent work on “control functionals” by Oates et al. (2014), we consider eval¬ 


uation of (p on two sets and at a computational cost (asymptotically) equivalent 
to evaluating ip on one such set. The first set is used to compute an arithmetic mean 
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based on a surrogate function : [0,1]'^ —)■ M. This surrogate function is itself estimated 
from the second set , in a preliminary step. In situations where ip^ can be made 
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to satisfy (i) f (pp^{u)du = f ip{u)du for all G N and (ii) V{(pjq) —)• 0 as iV —>■ cx), 
then the control fnnctional estimator Jcf is nnbiased (in an appropriate sense) and has 
asymptotically zero error relative to the standard QMC estimator. Oates et al.\ (2015) 


provides an explicit implementation of this strategy in the more general reprodncing kernel 


Hilbert space formnlation of QMC methodology (Dick et al, 2013) 


As a simple example, we note that for differentiable ip with snfhciently regnlar partial 
derivatives, a basic implementation prodnces a total variation V{(Pn) that vanishes at a rate 
Thus control functional QMC estimators are asymptotically superior to standard 
(R)QMC estimators under appropriate regularity conditions. Preliminary empirical results 
strongly support our theoretical analysis; an example is given in Fig. 

Given the gains in accuracy that are provided by QMC, it is surely a priority to establish 
complementary methodology that targets the rate constant governing the practical perfor¬ 
mance of these algorithms. Control functionals provide one (explicit) route to achieve this 
goal. The combination of control functionals with the Sequential QMC approach of Gerber 
and Ghopin should provide a highly effective approach to estimation. 


References 

Dick, J., Kuo, F. Y., Sloan, I. H. (2013) High-Dimensional Integration: The Quasi-Monte 
Garlo Way. Acta Numerica, 22:133-288. 

Hickernell, F. J., Lemieux, G., Owen, A. B. (2005) Gontrol Variates for Quasi-Monte Garlo. 
Statistical Science, 20(1):1-31. 

Oates, G. J., Girolami, M., Ghopin, N. (2014) Gontrol Functionals for Monte Garlo Integra¬ 
tion. CRiSM Working Paper Series, The University of Warwick, 14:22. 

Oates, G. J., Girolami, M. (2015) Variance Reduction for Quasi-Monte Garlo. Forthcoming. 

Owen, A. B. (1997) Scrambled net variance for integrals of smooth functions. Annals of 
Statistics 25 (4):1541-1562. 


2 








Figure 1: Convergence of control functional (CF) Quasi-Monte Carlo (QMC) estimators 
with respect to more standard alternatives. This example considers the (one-dimensional) 
test function (p{x) = sin(27ra;) -|- Ax. The usual Monte Carlo (MC) rate is and 

the usual QMC rate is (for any e > 0). The QMC-I-CF approach has a rate that 

is [Here we present the Randomised QMC (RQMC) case: For we used a 

scrambled Halton sequence of length N with a uniform random shift modulo one and for 
we used a uniform grid on the unit interval. Error bars denote ±1 standard deviation. 
For a fair comparison, the same number of function evaluations for (p was used in evaluating 
each estimator. It is striking how much the estimation error can be reduced by using control 
functionals.] 
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